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Recent work has shown that a simple chain of interacting spins can be used as a medium for 
high-fidelity quantum communication. We describe a scheme for quantum communication using a 
spin system that conserves z-spin, but otherwise is arbitrary. The sender and receiver are assumed 
to directly control several spins each, with the sender encoding the message state onto the larger 
state-space of her control spins. We show how to find the encoding that maximises the fidelity of 
communication, using a simple method based on the singular-value decomposition. Also, we show 
that this solution can be used to increase communication fidelity in a rather different circumstance: 
where no encoding of initial states is used, but where the sender and receiver control exactly two spins 
each and vary the interactions on those spins over time. The methods presented are computationally 
efficient, and numerical examples are given for systems having up to 300 spins. 



I. INTRODUCTION 

Quantum communication, the transfer of a quantum 
state from one place or object to another, is an important 
task in quantum information science [lj. The problem of 
communicating quantum information is profoundly dif- 
ferent to the classical case001- For example, quantum 
communication could not possibly be achieved by just 
measuring an unknown state in one place, and recon- 
structing in another. Rather, an entire system of source, 
target, and medium must evolve in a way that maintains 
quantum coherence. 

In this paper we consider an idealised system of inter- 
acting spin-1/2 objects, isolated from the environment. 
The aim is to use the system's natural evolution to com- 
municate a qubit state from one part of the system to 
another. The motivation is that such a system could 
be used as a simple "quantum wire" in future quantum 
information-processing devices. The most obvious config- 
uration to choose is a simple one-dimensional open-ended 
chain, with interactions between nearest-neighbour spins, 
in which case we want the chain's evolution to transfer a 
qubit state from one end to the other. The methods in 
this paper apply to this simple type of chain, and also to 
spin networks of arbitrary graph. 

A number of interesting proposals exist for quantum 
communication through spin chains. In Q, the ID 
Heisenberg chain was considered, with coupling strengths 
constant over the length of the chain and with time. The 
idea was to initialise all spins in the "down" state, except 
the first spin, which was given the state of the qubit to 
be sent. After the system was allowed to evolve, the spin 
at the far end of the chain would then contain the sent 
state, to some level of fidelity. Simulations were carried 
out for a range of chain lengths, and it was shown that 
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the fidelity was high only for very small chains. 

In 0, a ID spin chain with XY couplings was con- 
sidered. Here, the coupling strengths were constant over 
time, but were made to vary over the length of the chain 
in a specific way. Like Q, the first spin was initialised 
in the state to be sent, with all other spins initialised to 
"down" . It was shown that this scheme allows a perfect 
state transfer to the far spin site, for any length of chain. 

In 0, a scheme was presented for high-fidelity quan- 
tum communication over a ring of spins with nearest- 
neighbour Heisenberg couplings, using coupling strengths 
constant over the length of the ring and over time. The 
sender and receiver are located diametrically opposite to 
one another. The authors showed that excitations travel 
around the ring in a way that can be described using 
a concept from classical wave theory, the dispersion re- 
lation. Using this insight, they constructed a scheme 
where the sender, who controls several adjacent spins, 
constructs an initial state that is a Gaussian pulse hav- 
ing a particular group velocity chosen to minimise the 
broadening of the pulse over time. Using this state for 
the encoding of the |1) basis of the qubit message, and the 
all-down state as the encoding of |0), an arbitrary qubit 
can be sent with high fidelity over rings of any size, so 
long as the number of spins that the sender controls is at 
least the cube root of the total number. 

Motivated by the results in 0, we pose the following 
problem. Say we are given the Hamiltonian for a system 
of interacting spins, where the graph of the interactions 
is not necessarily a ring structure, but is completely ar- 
bitrary. Also, the strength and type of interaction along 
each graph edge is arbitrary (so long as total z spin is con- 
served). The sender Alice controls some given subset of 
the spins, and the receiver Bob controls some other given 
subset. How does Alice encode the qubit to be sent onto 
the spins she controls, in order to maximise the fidelity 
of communication? We know from [6j that the Gaussian 
pulse provides a near-optimal fidelity for the case of a 
Heisenberg ring (and is optimal in the limit of large ring 
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sizes). What about other shapes of spin network? Can 
we find a general solution? 

We provide a simple and efficient method for finding 
the maximum- fidelity encoding of the |1) message basis 
state, for a general z-spin-conserving spin system. (We 
assume that the encoding for the |0) basis state is fixed to 
the all-down state). So, unlike the schemes in 0, p|, and 
, which use systems with interactions that have specific 
strengths and conform to a specific graph, our scheme 
is designed to "make the most" of whatever arbitrary 
system is given to us. We give a numerical example of our 
method, for a system of 300 spins (where Alice and Bob 
each control 20 spins), showing a near-perfect average 
fidelity. 

We give a second scheme for increasing fidelity, that 
does not use encoding of initial states, but relies on Alice 
and Bob dynamically controlling the interactions on their 
control spins. Here, the number of control spins is fixed 
at two each for Alice and Bob. We give a straightforward 
method for deriving control functions, that give a fidelity 
(and communication time) equal to the values that would 
result if Alice and Bob had instead each controlled many 
more spins (with static interactions) and used the op- 
timal initial-state encoding scheme. This method has 
the combined benefits of being applicable to arbitrary z- 
spin-conserving spin-chains, yet having a fixed two-spin 
"interface" with Alice and Bob. We give numerical exam- 
ples, and plot the derived control functions, for a f 04-spin 
and a 29-spin system, showing a near-perfect fidelity in 
each case. 

In the remainder of this introductory section, we briefly 
describe the assumptions behind our schemes, and define 
our notation. Sec.[H]describes our method of deriving the 
optimal message encoding. Sec. IIIII describes our scheme 
for increasing fidelity via dynamic control. Concluding 
remarks are made in Sec. IIVI 

A. Assumptions and notation 

The solution presented in this paper relies on two main 
assumptions, which we now list. Firstly, the system 
Hamiltonian must commute with Z tot , which we define 
to be the z-component of the total spin operator 

-tot = ytot^ ^tot) = Y^a h (1) 

3 

where <jj is the vector of Pauli operators (a x ,a v ,a z ) act- 
ing on the j-th spin. The Pauli operators in the basis 
"down" | 4) and "up" | f) are 
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Secondly, the spin system must be initialised to the all- 
down state | I) ® • • • ® | J.), before the communication is 
carried out. Note that the schemes in H,[| and @ also 
make use of these two assumptions. It could be argued 



that the first condition is reasonable because it follows 
from rotational invariance. Of course, any external mag- 
netic field will destroy this invariance, and in particular 
any magnetic field which is not in the z direction will 
mean that the z-component of total spin is no longer 
conserved. The Heisenberg and XY interactions are ex- 
amples of interactions that conserve Z tot . The second 
constraint might be rather difficult to achieve in prac- 
tice. One possibility would be to apply a strong polar- 
ising magnetic field in the z direction, over the entire 
system, and let the system relax to its ground state. 

In the remainder of the paper, in place of the nota- 
tion | J.) and | t) for the eigenstates of a z , we will use 
the equivalent but more convenient notation |0) and 
A computational basis state of the system is defined to 
be one where each spin is in either a |0) state or a |I) 
state. Note that the computational basis states are all 
eigenstates of Z tot , and the eigenvalue has one of N + 1 
possible values, given by the number of |0)s minus the 
number of |l)s. So in a system of N spin-1/2 objects, 
we can break the state space into N + 1 subspaces of 
different well-defined z-component of total spin. We use 
n = 0, . . . , N, to denote these subspaces. is 
the eigenspace of Z tot that is spanned by the (^) com- 
putational basis states that have n qubits in the |1) state 
and the rest in the |0) state. 

Since the system Hamiltonian H commutes with Z tot , 
a state in U {n) will remain in TiS n > under the evolution of 
H . TiS ' is one-dimensional; it is spanned by the all-zero 
state |0)<g>. . . <8>|0). So this state is a stationary state of 
H. 



II. THE OPTIMAL ENCODING SCHEME 

Say that Alice wishes to send the qubit state a|0) + 
/3|f). In our scheme, she does so by preparing the state 
«|0) a+ P\^enc) A, where |0)^4 is the all-zero state on her 
spins, and \1enc)a is some state orthogonal to 10)^ (the 
"ENC stands for "encoded"). (Note that Alice doesn't 
necessarily know a and (3. She would presumably prepare 
the state by some unitary operation acting on her spins 
and some external spin containing the state a|0) + /3|I).) 
We assume that the entire spin chain is initialised to the 
all-zero state, so immediately after Alice prepares the 
abovementioned state on her spins, the state of the whole 
system is 

|*(0)) = (a\0) A + P\1 E nc)a) ® \0) A , (3) 

where A refers to all spins that Alice does not control. 
The whole spin system is allowed to evolve for a time T, 
giving the state |*(T)) = e" iHT |*(0)). Using the fact 
that |0) (g> • • • <8> |0) is a stationary state, \^(T)) can be 
written (up to some global phase) as 

|*(T))= Py/1-C B (T)\T,(T)) + 

|0} s H0) B +/VC^| 7 (T)) B ), (4) 



3 



for some nonnegative Cb(T), some normalised \j(T))b 
orthogonal to |0)s, and for some normalised \rj(T)) that 
is orthogonal to all states of the form \0)§ ® \v)b- 

We now show that Cb(T) can be used as a measure 
of success. Comparing Eqs. (0 and we see that 
Cb{0) — 0. If Cb(T) reaches 1 for some later T, a perfect- 
fidelity quantum communication has resulted. This is 
because Bob will then have the state a\0) b + (3\j) b on 
the qubits he controls, which can be "decoded" by a 
unitary operation into the state a\0) + of a sin- 
gle spin, since |0) B and \j) B are orthogonal. If Cb{T) 
is less than 1, the unitary decoding by Bob will leave 
him with a qubit state p that is generally different to 
the message state. That is, the measure of state fidelity 
F = (a|Q) + /5| 1> )+ p(a\0) + /3\1}) between the message 
a\0) +P\ 1) and p, will generally be less than one whenever 
Cb(T) < 1. However, the value of F is highly dependent 
on the message state — for example, if a — 1 then F = 1 
regardless of the value of Cb{T). 

Cb{T), on the other hand, is a message- independent 
measure of the fidelity of communication. Consider F, 
defined to be the state fidelity F averaged over all mes- 
sage states. For encodings \1enc) that belong to the 
Ti^ subspace, we have 
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r = § + 3VCs(T) + gC B (r), 



(5) 



which is a monotonic function of Cb(T) |(|. So, in this 
case maximising the average state fidelity is equivalent 
to maximising Cb{T). More generally, for \1enc) n °t 
m H {1 \ the expression in Eq. JSJ) provides a reasonably 
tight lower bound on F: 

l + iV^) + lcB(T)<F<l + iVCMf) + l 

2 1 



= - + -VCb(T). (6) 

The precise value of F will then depend on \r](T)) and 
the full specification of Bob's decoding unitary. 

A further argument for using Cb(T) as a measure of 
communication fidelity comes from considering the sys- 
tem's ability to transfer quantum entanglement from Al- 
ice to Bob. Suppose that Alice, instead of sending a mes- 
sage which is a pure quantum state a\0) + f3\l), sends a 
state which is maximally entangled with some additional 
spin that Alice possesses. (The additional spin does not 
interact when the system evolves). If the communica- 
tion is perfect, the result must be that Bob's decoded 
message becomes maximally entangled with Alice's ad- 
ditional spin. So more generally, when the communica- 
tion is not perfect, the amount of entanglement gener- 
ated between Alice and Bob would be a good measure of 
communication fidelity. In fact, the entanglement gener- 
ated, measured by the concurrence, is equal to y/C B {T) 
(a proof of this fact is outlined in Appendix A). This is 
independent of \r](T)), or the full specification of Bob's 
decoding unitary, or whether |1)_eatc belongs to 7i^\ 



To recap, when the Hamiltonian commutes with Z tot 
and the state is initialized to |0), the problem of achiev- 
ing a high communication fidelity can be boiled down to 
choosing an appropriate initial encoding \1enc) for the 
|1) qubit basis state. We seek a state \Ienc)a ®\®)a 
that has the property that it evolves to (or near to) a 
state of the form \0) g (g> \i)B, or in other words such that 
Cb{T) ~ 1 for some T. Alice's choice for the "encoding" 
of the |0) qubit basis state is fixed to the all-zero state. 
With perfect fidelity that basis state will evolve to the 
all-zero state on Bob's spins. (Note that in some cases it 
may be possible to increase fidelity further by allowing an 
encoding for |0) other than the all-zero state. We ignore 
such a possibility, in order to keep the method for finding 
the encoding simple and efficient. The simplification is 
used likewise in pfl.) 

We now show that the encoding \1enc) which max- 
imises Cb{T) for a given T can be found by performing 
the singular value decomposition on a modified version 
of the evolution matrix e~ lHT . 

Let A be the vector subspace of states of the form 
\^enc)a <8> |0)j[, such that a(0\^enc)a = 0. Similarly, 
let B be the vector subspace of states of the form |0)g ® 
\i)b, such that b(0|7)_b = 0. In other words, A reflects 
all the possible encodings that Alice could use for the |1) 
qubit basis state (regardless of the fidelity they would 
achieve). B is the set of states that we would like some 
state in A to evolve to; a state in A that evolves to one 
in B represents an encoding for |1) that gives Cb{T) = 1 
and thus a perfect average fidelity. 

Let P4 and P& be the projectors onto the subspaces 
A and B. Let U(T) = e~ lHT be the time-evolution op- 
erator. From Eqs. J2J and (@J, we can write 



Cb(T)= II P B U(T)\l 
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where |j • || denotes the ^-norm. This means that for a par- 
ticular total communication time T, choosing the optimal 
initial encoding for the |1) state is a matter of finding the 

normalised \ip) <E C 2 that maximises ||Pg[/(T)P4|'!/>)||. 
The maximum value is given by the largest singular value 
of U(T) = PbU(T)Pa, and the corresponding optimal 
|^i) is the first right-singular- vector of U(T) @- Recall, 
the SVD (singular value decomposition) of U (T) is 



U(T) = VSW^ 



(8) 
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where s\ > S2 > • • • > are the singular values, the 
orthonormal Wj are the right singular vectors, and the 
orthonormal ilj are the left singular vectors of U(T). Nu- 
merical packages such as Matlab have in-built routines 
for easily calculating the SVD. So, we have that Cb{T) 
has its maximum value, , when Alice chooses the initial 
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wi to encode |1). 
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Other parts of the decomposition could be useful as 
well. Say Alice wants to transmit two qubits simultane- 
ously to Bob. If she uses the all-down state to encode the 
1 00) basis state, then she should use the encodings w±, 
W2, and W3 for the other three basis states |01), |10), and 
1 11). Then, so long as S3 w 1, the two qubits would be 
simultaneously communicated with high fidelity. 

The vectors Wj and the values Sj are also the 
eigenvectors and square-root eigenvalues respectively of 
(P B U(T)P A yP B U(T)P A = P A UHT)P B U(T)P A . Now, 
Z tot commutes with P A U\T)P B U{T)P A because it com- 
mutes with each of P A , P Bl U(T), and W(T) separately. 
So the Wj will all have well-defined total Z spin (or can 
be chosen to, wherever ambiguities exist because of de- 
generacies in the Sj). This is important when it comes to 
calculating these solutions efficiently. Instead of perform- 
ing the full 2 N by 2^ matrix exponential and SVD, the 
calculation can be done separately for each of the smaller 
subspaces TL^ n \ starting each calculation with the (^)- 
by-(^) part of the Hamiltonian that acts on the TL^ 
subspace. 

Alice can't create a state with more than \A\ qubits 
in the "one" state, where \A\ is the number of spins she 
controls. So, in fact the calculation only needs to be done 
over the TiP^ , ... , H^ A ^ subspaces (in other words, the 
singular values corresponding to states in other subspaces 
will always be zero). 

In practice we have found that the optimal solution wi 
often belongs to the TiS^ subspace. (In particular, we 
calculated the optimal solution for a range of different 
values of T for various 8 and 9-spin systems, and found 
that only for a very small minority of the values of T, for 
each system, was the solution not in the Ti^ subspace ). 
A rudimentary argument for this can be made as follows. 
Looking for solutions in 'H^ n ' means optimising over Al- 
ice's ('J) degrees of freedom (of the space A n W (n) ), 
in order to make the final state land in or near a ( „ )- 
dimensional target space Br\H {n K This must be achieved 
despite the fact that the Hamiltonian is "trying" to move 
the state through a much larger (^) -dimensional space 
U [n) . Over the various values of n = 1, . . . , \A\, the di- 
mensionality of ACiTly 1 ' and BflTv- 71 ' as a fraction of the 
dimensionality of TL^ is largest when n = 1. In other 
words, when n = 1, the size of the target space, and 
amount of control available of the initial state, is largest 
as a fraction of the dimensionality of the entire subspace 

n (n) . 

So, in general we can restrict all the calculations to 
the V-dimensional subspace and there will still be 

a good chance that we will arrive exactly at the globally- 
optimal encoding wi . Ignoring solutions in the other sub- 
spaces will increase the efficiency of calculation consider- 
ably, especially for large chains. 

The evolution of a state in the TiS 1 ' subspace can also 
be interpreted as a continuous quantum walk of a par- 
ticle over a graph. ( For an introduction to quantum 
walks, see for example 9] and references therein). The 
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FIG. 1: The largest four sing ular values of U W (T), for a 
range of communications times T. Open-ended Heisenberg 
chain, iV=300, and Alice and Bob each control 20 sites. 



graph is simply the graph of interactions between spins in 
the Hamiltonian H, and the state \0}j <g) \l)j corresponds 
to the particle being at vertex j of that graph. So our 
methods for increasing communication fidelity are, equiv- 
alently, methods for guiding a quantum walk from one 
part of a graph to another. We point out this connec- 
tion because of the significant interest currently in using 
quantum walks for solving computational problems (see 
for example 0, ITol ITT| and references therein) . 

To demonstrate the use of the SVD optimal-encoding 
technique, we now consider a numerical example. Imag- 
ine that Alice and Bob are joined by a 300-site open- 
ended chain with nearest-neighbour couplings given by 
the antiferromagnetic Heisenberg interaction, with cou- 
pling strengths all equal to 1. That is, 

299 

# = I>;-^+i- (9) 
3=1 

Assume that Alice and Bob control the first and last 20 
spins respectively. 

In light of the earlier discussion, we restrict our optimi- 
sation to the subspace, and thus ignore all singular 
vectors in other subspaces. A Matlab program is used 
to carry out the following calculations. First, the 300 by 
300 matrix , defined to be the part of H that acts 
on Ti^\ is constructed. Then, the SVD of 

U {1) (T) = P mn(1) e^ H(1)T P Ann{1) (10) 

is calculated for a range of values of T. The optimal 
value for communication time is not known beforehand, 
so this repetition of the calculation for different values of 
T is needed in order to find a reasonable tradeoff between 
communication time and fidelity. 

The four largest singular values, 8%, .. . , S4, of U^(T) 
are plotted in Figure^ Over the range of T shown, si (T) 
has its maximum of 0.99999 at T = 75.75. So, this sys- 
tem can transmit a qubit with near-perfect fidelity, over 
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FIG. 2: The optimal encoded state wi (75.75), evolved for a 
sequence of times t. 



a time interval of 75.75. The graph shows that S2(75.75) 
and s 3 (75.75) are also very close to 1, so in fact two 
qubits could be transmitted simultaneously with high fi- 
delity in this example, using the encodings |0), wi (75.75), 
W2 (75.75) and wz (75.75) for the two-qubit basis states. 

Let's look at the actual optimally encoded states that 
are generated in this example. We visualise a state in 
TiP^ by plotting the square magnitude of the coefficients 
ipj , where ijjj is the coefficient of the basis state that has 
the j-th spin in the |1) state: 

n^3\^)= J2 (ii) 

• i v 

In Figure [3 we show the evolution of the state 
wi (75.75). That is, we set |^>(0)) = w x (75.75), and 

\ip(t)) = e _i - ff<1>t |?/'(0)), and plot the magnitudes \ipj\ 2 
for a sequence of equally-spaced times t. As is necessar- 
ily the case, the t = state has non-zero coefficients only 
on Alice's spins, j — 1 . . . 20. The state deforms itself 
into a Gaussian shape quite quickly. This is interest- 
ing in comparison with the results in Whilst Gaus- 
sian initial states were shown to optimise fidelity on a 
Heisenberg ring, the best initial states for an open-ended 
Heisenberg chain are ones that deform into a Gaussians. 
From the total communication time in this example, the 
group velocity of the pulse is roughly 3.95 (defining the 
distance between neighbouring spins to be 1). Thus, in 
the open-ended Heisenberg chain we have found the same 
phenomena that appeared in the Heisenberg ring in 
notably that the system has a preferred group velocity 
that gives a minimum dispersion and thus maximum fi- 
delity. This explains the fact that in Figure^the singular 
values drop for T greater than 75.75, and rise again to a 
near-maximum at 225 ~ 3 x 75.75: the high-fidelity 
communication for T 225 is also operating at the pre- 
ferred group velocity, but the wave packet is traversing 
the chain three times, after bouncing from each end. 

Curiously, the lower solutions W2 (75.75) and 103 (75.75) 
seem to evolve into a sum of two and three Gaussians re- 




FIG. 3: The states (a) w 2 (75.75) and (b) w 3 (75.75) evolved 
for time=12. 



spectively (see Figure |3J|. Animations of these evolutions 
are available online [12j. 

III. DYNAMIC CONTROL 

The scheme we presented in Sec. |n] utilises the evo- 
lution of a system having completely static interactions. 
The control that Alice and Bob have over the chain is 
only for an instant at the beginning and end of the pro- 
cedure, and so their only degrees of freedom for increasing 
the fidelity lie in the encoding they use. For very long 
chains, the number of control sites needed to give a high 
fidelity might become impractically large, as suggested 
by the results in 6j. In this section, we consider the ad- 
vantage that can be gained by allowing Alice and Bob to 
control their spins throughout the procedure, by modu- 
lating the strength of interactions on those spins. The 
advantages of the scheme are that the number of control 
spins are fixed at four, and that suitable functions for 
Alice and Bob to use to vary the interaction strengths 
are easily derived from a simple extension of the SVD 
approach already described. 

This type of control scheme is an example of a funda- 
mental problem in quantum information processing, that 
of determining how to use the limited physical control 
that one has of a quantum system, in such a way as to 
achieve the dynamics that are required. For the task at 
hand, our method provides a practical and efficient way 
of finding an appropriate dynamical control. 

A general schematic for the system is shown in Figure 
0] Alice is now in control of just two spins, labeled Ai 
and A%, and, likewise, Bob controls two spins B\ and £> 2 - 
All the other spins in the system are collectively denoted 
C. The graph of the interactions connecting Alice and 




FIG. 4: The general setup, whereby Alice and Bob modu- 
late a total of six parameters of the Hamiltonian in order to 
increase communication fidelity. 
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FIG. 5: This is a modified version of Figure Calculating 
the optimal initial encoded state on this system will help Al- 
ice and Bob derive suitable control functions for their actual 
physical setup in Figure @. 



Bob's spins can be arbitrary, except that A\ must directly 
couple only to A2, and B\ must directly couple only to 
£>2 • Like the previous scheme, we require that the system 
Hamiltonian H commutes with Z tot , and that all spins 
are initialised to |0) before the procedure starts. So again 
the problem is that of finding a way of sending the |1) 
basis state with high fidelity 

The protocol works as follows. At time t = 0, Al- 
ice transfers the qubit state she wishes to send onto the 
spin A\. Then she varies the coupling strength between 
A\ and A2 according to some function which we de- 
note Ja(£) ; and varies the z magnetic field on A\ and 
A2 according to functions BAi(t) and BA2(t). At the 
same time, Bob varies his coupling strength and z mag- 
netic fields according to the functions Jp,(t), BBi(t) and 
Bp2(t)- This process is continued over some time inter- 
val < t < T. At t — T Bob's spin B Y will contain the 
sent qubit state, to a level of fidelity that depends on the 
six control functions. 

How do we choose the control functions in a way that 
gives us a high average communication fidelity? The trick 
is to imagine a modified version of the system, where a 
number of "phantom" spins have been added to both 
Alice and Bob's set of control spins, but all couplings are 
now fixed (see Figure [5J. 

The SVD method is applied to this modified system, 
to find the optimal initial state on Alice's extended set of 
control spins. The evolution of the encoded state through 
the modified system is then simulated on a (classical) 
computer, and the results of the simulation are used to 
derive appropriate control functions for the actual phys- 
ical system, using a method which we describe below. 
Since, over the bulk of the physical system, the initial 
state and interactions are identical to those at the cor- 
responding regions of the modified system, the problem 
reduces to finding control functions which make the state 
on A2 and B2 in the physical system evolve in the same 
way as those corresponding spins in the modified chain. 
When that is achieved, the state on the bulk of the phys- 
ical chain will evolve in the same way as in the modified 
chain, which means that we can communicate a qubit 
with the same fidelity as for the optimally encoded state 
in the modified system. 

For the sake of clarity, we describe the method in detail 
for a less general configuration, the ID XY chain. The 



derivation is simpler in this case because, as we shall see, 
the magnetic control is not needed. The system Hamil- 
tonian is given by 

JV-2 

ff(*) = E J i[°i°i+i + ^°3+i] 

i=2 

+J A (t) (a^ + af4) 
+J B (t)(<T x N _ 1 a* N + *l_ 1 a%). (12) 

That is, there are some arbitrary fixed Jj that specify 
the strengths of the XY couplings over the bulk of the 
chain. The strengths at the first and last links can be 
varied over time by Alice and Bob. 

We write down a Hamiltonian if of a modified system, 
where we have extended the length of the chain in both 
directions by adding Np phantom spins to both Alice 
and Bob's sides. The new coupling strengths are chosen 
to be 1, and the two strengths that were time- varying in 
the original system are now also fixed at one. So that 
we can use the same numbering system for the spins as 
in Eq. <|12|l . we let the indices of the spins range into the 
negative in the modified system, running from 1 — Np to 
N + Np. The modified Hamiltonian is written simply as 

N+Np-l 
] = 1-N P 

where we have extended the definition of Jj so that it 
equals 1 for j = (1 — Np), . . . , 1 and for j — (N — 
1),...,(A + A P -1). 

Next, using the SVD method, we find the best encoded 
initial state on the set of spins from index (1 — Np) to 1, 
while assuming that the target set of spins ranges from N 
to (N + Np — 1). That is, we are imagining a "modified 
Alice" that controls the first (Np + 1) spins and a "mod- 
ified Bob" that controls the last (Np + 1) spins of this 
extended chain. Recall that the SVD method depends on 
a choice of total communication time T. As in the exam- 
ple earlier, we may wish to search over a range of values 
of T to find the most suitable value. It is important for 
the procedure at hand that we restrict ourselves to states 
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in the liS 1 ^ subspace (whereas before this restriction was 
just a way of making the solution much faster to com- 
pute). So, we calculate wi(T), the first right-singular- 
vector of PBr\nw e ~ 1 ^ {1)T PAnnw > where is the part 
of H that acts on the subspace. 

Then, we need to be able to calculate the evolution of 
the state Wi, over a range of times t from to T. Let 



|V(0))=«Ti(T), and |^(t)) = 



|?A(0)). As earlier, 



the evolving state is a series of complex coefficients i/jj (t) , 
where j is the index to a spin site, ranging from (1 — Np) 
to (N + Np-1). We need to know tpi(t) and tp N (t) for 
every value of i that we wish to calculate J Ait) and J sit) 
for. 

Similarly we use 4>j{t), j — 1, . . . , N to denote the evo- 
lution of the |1) qubit state over the original physical 
chain. Recall that in this scheme, Alice places the qubit 
state to be sent, unencoded, onto spin number 1, after 
all other spins have been initialised to zero. So, initially 
we have <fo(0) = 5jA- The functions (j>j(t) depend on the 
control functions Ja(£) and Js(t) (whereas the ipj{t) do 
not). 

The aim is to chose control functions Ja(^) and Js(t) 
in such a way as to force 4>j it) = ipj it) , for all the spins 
in the range j = 2, . . . , N — 1, and for all t in the interval 
[0, T]. That is, we know the way the optimal encoded 
state evolves over the modified chain, and we want to 
make the |1) state in the physical system evolve in exactly 
the same way, over all spins except 1 and N. In this way, 
the physical system will carry a qubit across it's length 
with the same fidelity as the encoded modified system 
does. 

The interactions on the spins from site 3 to site N — 2 
are the same in the physical chain as in the modified 
chain. So, the differential equations for the are 
the same as those for the 4>jit), f° r j — 3, ... ,N — 2. 
Specifically, 



dt 
dt 



= —2i [Jj_i<ipj_i{t) + Jjipj + i(t)] and 

= -2% [Jj-x^-iit) + Jrfj+iit)} , (14) 



for j = 3, . . . , N — 2. Also, the initial conditions are the 
same between the ipj and the <f)j, for j = 2, . . . , (N — 1): 
^(0) = ^(0) = 0. 

It follows that if we can use our control functions to 



force 



dipN-i(t) _ d0jy_i(t) 



the 



time range t = 0, . . . ,T, then we will have = 4>jit) 

for all t in that time range, and for all j = 2, . . . , N — 1, 
as desired. 
Now, 



#2 it) 

dt 



-2i[Mt) + J2Mt)} 



(15) 



and 



So, assuming that at time t <j>j{t) = i/>j(t) for j 
2, . . . , N - 1, then = ^1 by setting 



Mt) 



and 



i(t) _ rfiAiv-i(t) 



dt 



dt 



by setting 



Jsit) 



ip N (t) 
it)' 



>N 



(17) 



(18) 



Thus, the practical task of numerically calculating 
Ja(£) and Je(i) involves simulating the evolution of both 
the <pj and ipj states on the original and modified systems 
respectively, over the time interval [0, T], and evaluating 
Eqs. ((n|| and iPfl. 

The functions J a it) and Jb it) must of course be real- 
valued, for the Hamiltonian to be Hermitian. Equations 
(fTTjl and {TBI will indeed be real for the XY chain. The 
expressions for the -$f are all given by a purely imagi- 
nary linear combination of the nearest-neighbour values 
4>j-i and 4>j+i- Then, considering the initial conditions, 
</>j(0) = <5i j, it's clear that the <f)jit) are real for odd j 
and imaginary for even j, for all values of t. The val- 
ues ipj(t) also have this property of alternating real and 
imaginary values. Again, the time derivatives of ipjit) 
are purely imaginary linear combinations of the values 
tjjj-i(t) and ipj+i{t). Thus, by performing the change of 
variables 



m) = 



ii)j{t) 



if j is odd, and 
if j is even, 



(19) 



the differential equations for ipjit) will all have real 
coefficients. So, the entries of the evolution matrix 
e -iH [ t mus t k e rea ^ after that change of variables. 

Thus, so must be the entries of PjsnHM & 1 ^-FUnftC 1 ) • 
So, w(T), which is the right-singular vector of 

^BrrHW e ~ lH 1 T -^UnW( 1 ) > wm a l so have all real coeffi- 
cients with respect to the changed variables. Changing 
variables back, the initial encoded state V'j(O) = WjiT) 
will thus have the property of having real values for odd 
j and imaginary values for even j, and so will ipj it) for all 
t. So, Eqs. (|17fl and l|18|) will be real-valued as required. 

Eqs. (T7J| and {TSJ will never be infinite. In fact, | Ja(£)| 
and |Js(t)| will be at most 1. This is a simple conse- 
quence of conservation of probability. Since ipj it) = <pj it) 
over the bulk of the chain (J = 2, . . . , N — 1) and Alice 
and Bob's sides only interact via the bulk of the chain 
for both the physical and modified systems, we have that 



J2 |^| 2 = l0i| 2 ,and 

N+N P 



E IVyl 3 = l* 

j=N 



N 



(20) 



(21) 



dMt) 
dt 



from which it follows that |'0i(t)| < |<^i(i)| and 1^(^)1 < 
-2i[J A it)Mt) + J2Mt)}- (16) \<p N (t)\. So, |J A (t)| < 1 and \ J B (t)\ < 1, if they are 
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defined. If Jyi(i) (or Js(i)) is undefined (0/0), it means 
that the requirement of = ( res P ec tively 

^"df 1 ^ = d '^ N dt^ > ) ^ S satisfied regardless of the value 
of J^(i) (respectively Jg(t)) for that i, in which case the 
value of the control function can be chosen arbitrarily at 
that time. 

We now plot the derived control functions Ja(^) and 
Js{t) for two simple example XY chain systems. We 
used numerical integration in these examples, in calcu- 
lating the evolution of the 4>j (t) due to the time- varying 
Hamiltonian. We divided the total evolution into a num- 
ber of discrete time steps, where the approximation was 
made that the Hamiltonian remains constant throughout 
each step. The value of Ja(0 and J sit) for a step was 
calculated from the state of the system at the previous 
step. We used 2000 time steps, which gave a final fidelity 
in the physical chain within two significant figures of the 
correct value given by the evolution of the static modified 
system. 

The first example is a chain 104 spins long (ie. 100 
non-controlled spins, plus the four control spins), with 
all the non-controlled coupling strengths set to the same 
value, 1. The control functions were derived by using 
a modified chain 144 spins long (that is, 20 phantom 
spins added to each side) with all coupling strengths set 
to 1, and the total communication time T chosen to be 
36. Figure © shows that the resulting J^(i) and J sit) 
are quite simple and well behaved. The fidelity measure 
Cb{T) is 1.0, to 6 decimal places. This can be compared 
with the fidelity in the same 104-spin system but without 
the time-dependent control, that is with J a (t) and Jb (t) 
fixed at 1: over the time interval < t < 1000, the value 
of C B {t) is at most 0.2809. 

The second example is an XY chain 29 spins long, 
but where the non-controlled coupling strengths are ran- 
domly sampled uniformly from the interval [0.95,1.05]. 
This is as if the chain has been manufactured with ran- 



dom imperfections in the coupling strengths, but these 
coupling strengths have been somehow measured after 
the manufacturing process and are known to Alice and 
Bob. (A shorter chain was chosen in this example, com- 
pared with the previous example, in order that a near- 
perfect fidelity would still result. We have observed that 
when random couplings are used, the achievable fidelity 
will decrease as a function of the chain length). We de- 
rived control functions using a modified system with 25 
phantom spins added to each side, where the new arbi- 
trary coupling strengths are set to 1. The communica- 
tion time T was chosen to be 19.5. Figure (JJJ) shows 
control functions which are a little more complicated in 
this case, but still rather smooth. The fidelity measure is 
Cb{T) = 0.99625. In comparison, in the non-controlled 
version of this system, with Ja^) = </b(£) = 1, the 
value of Cs(i) does not exceed 0.496 over the interval 
< t < 1000. Animations of both examples are available 
online [12j. 

What about a system that is not simply an XY chain, 
but any configuration conforming to Figure Q and con- 
serving Z tot l Then, the ideas and methods are almost 
the same, but with the added complication that we need 
control the z magnetic fields on Alice and Bob's qubits as 
well as controlling Ja(£) and Js(i). The are two reasons 
why the magnetic control is needed, which we explain 
for Alice's side. First, the simple phase relation between 
ipi(t) and (f>i(t) that we saw in the XY chain does not 
occur in general. So, BAi{t) is chosen simply to keep 
<j>\ (t) in constant relative phase to the ipi(t). Second, the 
type of interaction that the J a (t) is modulating may con- 
tain it's own magnetic-field-like interactions (that is, non- 
equal diagonal elements in the Hamiltonian) that need to 
be cancelled by Bai(£) and Bai (t). General expressions 
for BAi(t) and BA2(f) are straightforward to derive, but 
not particularly illuminating, so will not be given here. 
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IV. CONCLUSION 

We have considered the problem of communicating a 
quantum state over an arbitrary Z tot -conserving spin sys- 
tem. Our first scheme used a static system Hamiltonian, 
and utilised the fact that the sender and receiver con- 
trol several spins each, to increase fidelity by performing 
state encoding. We showed that choosing the optimal 
state encoding is a simple matter of performing a SVD 
on a modified evolution matrix. 

We have also shown that if the sender and receiver 
have control of just two spins each, but can vary the in- 
teractions on these four spins over time, then they can 
achieve a fidelity that is equal to if they each controlled 
many more spins on a static system and used the optimal 
state encoding. We have given a practical method of de- 
riving suitable control functions. The advantage of this 
scheme is the "fixed interface" that Alice and Bob have 
with the chain. That is, if the chain is altered, the only 
change that Alice and Bob need make is to their con- 
trol functions, rather than to the number of spins they 
control. 

It should be noted that the systems we have considered 
are idealised to a high degree. In particular, we haven't 
considered the effects of external noise, or the effect of 
having a Hamiltonian that only approximately commutes 
with Z tot , or the case where Alice and Bob have only 
an approximate knowledge of the system Hamiltonian. 
These issues will be the subject of future work by the 
author. 
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APPENDIX A 

Here we outline a proof of the claim in Sec. [H] regard- 
ing the connection between Cb (T) and the system's abil- 
ity to transmit entanglement from Alice to Bob. This 
connection helps establish Cb(T) as a good measure of 
communication fidelity. 

Proposition. Suppose that Alice sends a state which is 
maximally entangled with some additional spin that Al- 
ice possesses. (After the maximally entangled state is 
created, the additional spin is assumed to not interact 
during the remainder of the communication procedure). 
Then, after the communication procedure of Sec. El * s 
carried out, the entanglement (measured by concurrence) 
between Alice's additional spin and Bob's decoded mes- 
sage, equals ^/Cb{T). 

Proof: Note that it doesn't matter which maximally- 
entangled state is used — all such states are equivalent up 
to a local unitary on the additional spin, and such a local 



unitary could not possibly affect the way entanglement 
is transferred through the system. 

Let the additional spin "+" and the spin "M" con- 
taining the message have the maximally entangled state 
>+|1)m) 



|*(0)) 



) A 0) 



1J1 



+ \±ENC/A 



10) 



(Al) 



(|0) + |0)a/ + 1 1> + 1 1)a/)- Thus, after Alice performs her 
encoding, the entire state is: 

1 

After the system evolves for time T, the state becomes 



\0) S \0) B + \l) + Ml-C B (T)\r)(T)) 



(A2) 



Then Bob performs a decoding unitary, denoted t/dec, on 
the spins he controls, t/dcc is defined to act as follows: 
U dec \0) B = \0)b and U dec \^(T)) B = |0...01} B , where 
|0...01)s is the |1) state on spin N and the all-zero 
state on Bob's other spins. After Bob's decoding, the 
joint state of Alice's additional spin and Bob's decoded 
spin is: 

P+/N = tr w (C/ dec tr s (|$(T))($(T)|)Z7] eG ) 

= -C B (T))\1){1\ ® p + C B (T)|ll><ll| 



-|00)(00| + VCb(T)(|00)<11| 

-IhXoodI, 



(A3) 



where p = tr T7 ^(J7 d octr i3 (|?7(r))(r/(T)|)C/] cc ), and where 

tr(.)(-) is the partial trace performed over the spins indi- 
cated. 

Concurrence is a measure of entanglement between two 
qubits jl3l | . The value of concurrence for a density matrix 
P+/N is equal to 



e (p+/n) 



c{0, 



A 4 }, (A4) 



where the Aj-s are the eigenvalues, in nonincrcasing or- 
der, of the matrix p + / N (a y ® a v )p* + ^ N (a y ® a v ), where 
* represents complex conjugation in the computational 
basis. It can be shown that 



Ai = ~ ( v / /5iiC B (T) + l- i 5 11 + y/C B (T) 
A 2 = \ (^JpuCb{T) + 1- pu - VCb{T) 



A 3 
A 4 



(A5) 



where pn = (0\p\Q). Thus, using the fact that Cb(T) and 
pn each lie in the interval [0, 1], we have 



E(p +/N ) = y/C B (T), 



as required. 



(A6) 
□ 
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